Soliton stability in a generalized sine-Gordon potential 
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Abstract 

We study stability of a generalized sine-Gordon model with two coupled scalar 
fields in two dimensions. Topological soliton solutions are found from the first-order 
equations that solve the equations of motion. The perturbation equations can be cast in 
terms of a Schrodinger-like operators for fluctuations and their spectra are calculated. 

PACS numbers: 0.350.-Z, ll.lO.Ef 
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1 Introduction 

It is well known that in field theories when a discrete symmetry is broken domain walls arise. 
Domain walls have been observed in condensed matter, for example, in liquid crystals. In 
the cosmological context domain walls could appear in phase transitions in the early universe 
and have some important consequences (Vilenkin and Shellard, 1994). 

In the domain walls context, there exist classical static configurations with finite mini- 
mum localized energy, see for instance (Bogomol'nyi, 1976; Prasad, 1975). Several authors 
have been interested in coupled scalar fields systems due to their important physical prop- 
erties. For example, Peter showed in (Peter, 1996) that surface current-carrying domain 
wall arises when a bosonic charge carrier is coupled to the Higgs field forming the wall. In 
(Bazeia, et al., 1997; Riazi, et al., 2001) was studied linear stability of soliton solutions for 
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a class of systems of coupled scalar self-interacting fields following the standard approach 
of classical stability. 

Besides, as it was extensively shown, Witten's supersymmetric quantum mechanics (Wit- 
ten, 1981) is the most recent way to study solvable as well as perturbative problems (Witten, 
1981; Lahiri et al., 1990; Cooper et al., 1995; Junker, 1996; Cooper et al, 2001). 

Some systems of coupled scalar fields present soliton solutions and their linear stability 
have been addressed by means of the SUSY QM formalism in a 2 x 2 superpotential realiza- 
tion (Bazeia et al., 1995; Bazeia and Santos, 1996; de Lima Rodrigues et al., 1998; Dias at 
al., 2002). The superpotential satisfies the Riccati equation associated to the perturbation 
Hessian and establishes the stability of the system. This fact is because of the correspond- 
ing supersymmetric operators factorize the perturbation equation and automatically ensure 
non-negative equal perturbation frequencies (Bazeia and Santos, 1996). 

The sine-Gordon system has been applied to a wide class of physical problems like prop- 
agation of crystal dislocations, two-dimensional models of elementary particles, propagation 
of splay waves in membranes, Bloch wall motion in magnetic crystals and magnetic flux in 
Josephson lines (Rajaraman, 1982). It is well known that this system in (1-1-1) dimensions 
has classical soliton solutions and their non-dissipative properties could be explained like a 
finely-tuned balance between self-interactions and dispersion. 

A generalization of sine-Gordon system with two coupled real scalar fields showed an 
important and rich behavior (Riazi et al., 2002). The potential of this system consists of 
a product between a trigonometric and polynomial functions of the fields. Depending on 
the rest energies and the boundary conditions, the spectrum of solitons could be stable, 
unstable or meta-stable. The former classical stability analysis was established by means of 
numerical analysis. 

Another generalization of coupled sine-Gordon model has been given as an example of 
continuously degenerate soliton (Shifman and Voloshin, 1998). In contrast to the mod- 
els mention above, this model involves a highly coupled self-interacting fields with non- 
polynomial form. This model has richer structure and dynamics and deserves further anal- 
ysis. 

In this paper we advocate to the study of linear stability approximation of the gener- 
alized sine-Gordon model proposed in (Shifman and Voloshin, 1998). In section 2, we give 
the model consisting of two couple real scalar fields, we present the first order differential 
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equations that minimize the energy for the fields and find a particular solution of them. 
In section 3, we show that the stability equations can be analyzed in terms of SUSY QM 
formalism and reduce the problem of stability to solve a Riccati equation associated to the 
perturbation Hessian. In section 4, we find the spectra of the fluctuation operator and 
explictly show the stability of the soliton solution. Finally, we give the conclusions of this 
work. 



2 The Model 



In this paper we consider the generalization of the sine-Gordon model for two scalar inter- 
acting fields given by the following Lagrangian 



"^'""^2 icos2$(l + asinX)2 - icos2X(l + asin$)2, (1) 



where a is a dimensional parameter, and all other dimensional parameters are set equal to 
unity. For < $, X < 27r, the potential in ([T]) has three minima <^ ~ X — tt/2; $ = 7r/2, 
X ~ 37r/2 and <f> X = 37r/2, one maximum &i ^ ~ Q, X — tt /2, and three saddle points 
at $ = it/2, X = tt; $ = 7r/2, X = and $ = 37r/2, X = 7r/2. 
The equations of motion for the model ([!]) are the usual ones: 



□$-(-— -V = 0, 

which become for a static configurations 

<i>" = -cos$sin$(l + asinX)^ - 
X" = -cosXsinX(l + asin$)2 



acos^ X{1 + asinX) cos<I> 
- acos^ $(1 + asin(f>) cosX, 



(2) 

(3) 
(4) 



where primes means derivatives with respect to space variable. 
The form of the energy of the system can be written as 



— -W^ 

az 



dX 
dz 



W 



X 



dz 



" ^WMz),X{z)]dz 



(5) 



where X(z)] is the corresponding superpotential of ([T]), which turns out to be 

= - sin* - sinX - Q!(sin<I>)(sinX). (6) 

In (Shifman and Voloshin, 1998) this superpotential was referred as a generalization of the 
sine-Gordon model. It is periodic in both $ and X; for a = it describes two decoupled 
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fields, representing each of them a supergeneralization of the sine-Gordon model. If a 7^ 
the fields <& and X start interacting with each other. Inside the periodicity domain < 
$, X < 27r, —W has one maximum at $ = X = 7r/2, one minimum at $ = X = 37r/2 and 
two saddle points at <& — 7r/2, X = 37r/2 and $ = 37r/2, 7r/2 at least for small values of a. 
The lower bound for the energy is achieved if $ and X satisfy 

$' = - cos<l)(l + asinX) 
X' = -cosX(l + asin$). (7) 



For the case X = 7r/2, we have 



— = -cos$(l + (8) 



whose solution is 



/g2g2:(l+Q) _ -z(l+a) - 

* = -tan-M ). (9) 



Other possible solution of equations ((?]) is obtained for X = 3tt/2 and the solution is obtained 
from the former equation by substituting a by —a. Interchanging the fields X and $ in the 
last equations we get the solution for <f> — tt/2 or $ = 37r/2. 

We attempted to find the general solutions of the coupled equations ([7]) by the trial orbit 
method of Rajaraman (Rajaraman, 1977), however, because of the difficulty of the system 
we were unable to find them. 



3 Stability Equations and SUSY QM 

We are interested in determining the classical stability of this system under small fluctuations 
around a static configuration. In order to investigate the linear stability of the interacting 
fields we proceed in the usual way by considering small perturbations around the static 
scalar fields 

= <^>{z)+rj{z,t) (10) 

x{z,t) - x{z)+az,t). (11) 

The stability equations can be written in a Schrodinger-like equation 

Sl^n = tOn^n (12) 
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where n = 0, 1, 2... The differential operator Si is given by 



c _ I cU^ ^ " d^dx I — T -L \7 ^^^^\ 



and the two components wave functions are 



|*=*(^),X=X(2) 




*n= , (14) 



where we have expanded the fluctuations o.(z,t) and P(z,t) in terms of normal modes 

r,{z,t) = J2''nVn{z)e'^"* (15) 

n 

az,t) = X^6„C„(z)e*'^"*. (16) 

n 

Notice that in the case when the differential operator Si is diagonal the perturbation fields 
could be expanded in terms of different frequencies. 

The SUSY QM approach to linear stability consists in realizing a 2 x 2-matrix super- 
potential, which is obtained by solving the Riccati equation associated to the perturbation 
Hessian Vph 

W^+W' = VpH. (17) 

The existence of W that satisfies this equation ensures the existence of the first order self- 
adjoint differential operators 

X>± = ±1-^ + W(z) (18) 

that factorize the operator Si = I>+I>~. This fact implies the stability for equal fluctuation 
frequencies, since < |X>-*„|2 = (X>-*„)t(D-'J'„) = (D+X>-) = {Si) = uj^. 
For our case the matrix elements of VpH are given by 

(Vph)!! = -(cos^$-sin24>)(l + asinX)2-|-Q!2cos^Xcos^$ (19) 

— acos^X(l + Q;sin$)sin$ 

(Vph)22 = -(cos2X-sin2X)(l+asin$)2+a2cos2$cos2X 

— acos^ $(1 + asinX) sinX 

(Vph)i2 = (Vpij)2i = -2acos$sin$(l + asinX)cosX (20) 

— 2Q!C0sXsinX(l + asin$) cos$. 
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The solution of Riccati equation (fTT]) for configurations satisfying equations ([7|) is 

(1 + a sinX) sin $ — acos$cosX 
-a cos $ cos X (1 + asin<i>) sinX 

For the sector X — tt/2, the fluctuation potential term becomes 

-(l + a)2(cos2$~sin^$) 

(1 + Q:sin<i>)2 - Qcos^ $(1 + a) 

so, the corresponding superpotential is 

(1 + a) sin$ 





(21) 



(22) 



(23) 



We point out the existence of another self-adjoint and non-negative second-order differ- 
ential operator S'^ = V^D^ which plays the role of the supcrsymmetric partner operator of 
Si in SUSY QM. The operators Si and S[ have the same energy spectrum except for the 
ground state. 



4 Spectrum of the Second Order Fluctuation Operator 

The study of stability for the general case is very difficult. However, in order to have 
analytical results in the case oi X — tt/2 (the results we are going to obtain are auto- 
matically true for $ = 7r/2), wc take the particular case of c = 1 in equation ([9]) i. e. 
tan<i> = sinhz(l + a). Since the differential operator Si is diagonal we could have different 
fluctuation frequencies that can be determined from the perturbation equations 

- ^ - (1 + a) (2sech2z(l + a) - l) r,„ = c^^^^ (24) 

and 

--^ + (l + - 2atanhz(l + a) ~ a(l + 2a)sech^z(l + a)) C„ = ul^n- (25) 

Performing the variable change y = z(l + a), equation transforms to the Rosen-Morse 
equation (Morse and Feshbach, 1953). We find that the fluctuation frequencies are 

u;l = {l + ar{l-il-nrf. (26) 
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However the bound states exist only for n < 1 (Morse and Feshbach, 1953). Thus, the 
ground state 770 = (1 + a)sechz(l + a) with eigenvalue cjo — is stable. 

By means of the same variable change the equation can be cast as a Rosen-Morse 
equation whose eigenvalues are 

Lol = l + a^-{l + af 
4a 

(1 + ")^ (5(T^^-(" + l/2))' 
which are the frequencies for possible bound states. However, for this case we have no bound 
states because n must be less than zero for both a > and a < (Morse and Feshbach, 
1953). 

This means that the solution configurations are stable under small perturbations around 
X = tt/2, tan<i> ~ — sinhz(l + a) ( the same is true for $ = tt/2, tanX = — sinhz(l + a)). 

Conclusions 

We have applied the SUSY QM formalism to study the linear stability of the Shifman gener- 
alization of the sine-Gordon model. We have shown that stability for soliton configurations 
is ensured by solving the Riccati equation for the 2x2 superpotential associated to the 
non-diagonal perturbation Hessian. The spectrum of the second order fiuctuation operator 
for the general case is very difficult to find it. Thus, we have got the fluctuation spectrum 
for the particular case X ^ tt/2 (or <f> = 7r/2), and we have found analytical solutions and 
explicitly found that the system is stable. We notice that equations (|24|) and ([25)) can be 
reduced to a Rosen-Morse equation. On another hand, the Rosen-Morse equation have been 
studied from the shape invariance approach of SUSY QM (Dutt et al., 1988). Thus, each 
one of the equations and ([25)1 has a scalar superpotential. Therefore, we have given a 
complete treatment of linear stability for the generalized sine-Gordon superpotential from 
the point of view of SUSY QM. 
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3a + 1 
2(1 + a) 



(n + 1/2) 



(27) 
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